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By Stanley U. Benscoter and Richard H. MacNeal 

SUMMARY 


An application is descrlhed of the well-known analogy bet’vreen elec- 
trical and mechanical systems to the calculation of stresses and deflec- 
tions of beams. The object of the present paper is to give an explanation 
of the analogies in an elementary manner vdilch will enable a structural 
engineer to londerstand the process of designing ’the electrical circuits. 
The analogies which are discussed ^e those that are now being used in 
the Cal-Tech analog computer. Analogies are given for beams in bending 
and torsion with static loads and in vibrational motion. 


INTRODUCTION 


In recent years a considerable amount of progress has been made in 
the development of applications of electrical-circuit analogies to prob- 
lems in structural analysis. Many solutions of such problems have been 
obtained on the Cal-Tech analog computer, A number of papers have also 
been written in which the electrical analogies for various problems have 
been presented. These papers have been written mainly by electrical 
engineers and have assumed the reader to be reasonably familiar with 
circuit theory. Consequently structural engineers have found some dif- 
ficulty in appreciating and evaluating this new method of analysis. In 
the present paper the design of the electrical circuits will be explained 
in complete detail without assuming any knowledge of circuit theory on 
the part of the reader. It is hoped that this paper will assist structural 
engineers in reaching a better understanding of the practical utility of 
electrical analogies in structural analysis. 

In order to develop an electrical analogy for the action of a beam 
under load it is necessary to replace the differential equations which 
govern the stresses and deflections by difference equations. The rela- 
tions between voltages and currents in an electrical circuit may also be 
expressed by difference equations. An analogous electrical circuit for 
a beam is one in which the voltages and currents are related by difference 
equations which have a form Identical to that of the structiiral equations. 
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including the boundary conditions. When such a circuit is designed and 
constructed the physical circuit, with its components, is called an 
analog computer. It might also be called an electrical simulator or 
even an electrical model. 

This paper was prepared at the California Institute of Technology 
under the sponsorship and with the financial assistance of the National 
Advisoiy Committee for Aeronautics. 


SYMBOLS 
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shear-resistant area 
capacity of a condenser 

Young's modulus j also voltage with respect to ground 
shearing modulus of elasticity 

moment of inertia of cross section of a beamj also current 
rotary inertia of a segment of a beam 

moment of inertia of mass about axis of beam 

current 

torsion constant 

station number along a rlbj also number of a nodal point 
station number along a beamj also number of a nodal point 
Inductance of an induction coll 
bending moment 

massj also bay number along a rib and number of a nodal point 
when used as a subscript 

bay number along a beamj also number of a nodal point 
concentrated vertical load 
concentrated torque load 
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Vsn 

Vin 
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X 
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resistance 

turns ratio of a tr^sformer 

section torque 

time 

section shear 

vertical component of internal forces in nth hay 
indeterminate part of section shear 
vertical deflection 
spanwise coordinate 

angle of rotation of a normal to elastic axis 

shearing strain 

differencing operators 

slope of elastic axis 

hay length 

Poisson’s ratio 

rotation (angle of twist) 

quantities in upper coil 

GENERAL DESIGN PRINCIPLES 


From a mathematical viewpoint the theory of circuit analogies has 
not heen fully developed. Questions of existence and uniqueness remain 
yet to he answered. From an engineering standpoint the procedure for 
designing a circuit which is analogous to a given structure is at present 
dependent to some extent upon the ingenuity of the designer. It would he 
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desirable, to have available a set of rules defining a standardized 
design procedure which could always be depended upon to give a correct 
analogous circuit for any given differential equation. However^ no such 
rules have as yet been developed. 

Several principles have been found to be useful in a large number 
of circuit designs and these will be stated briefly. Exact limitations 
upon the usefulness of the principles are not known. The first princi- 
ple is that the equations governing the action of the structure must be 
expressed as difference equations with respect to the space coordinates. 
If the stresses or deflections are governed by differential equations 
these equations must be replaced by difference equations. 

The two quantities which are available in an electrical circuit are 
voltages and currents. The two types of quantities which occur in struc- 
tural analysis are force quantities and displacement quantities. Two 
general methods can be used in designing analogous circuits. In one 
method the voltage drops across various elements of the circuit corre- 
spond to forces and the currents in the elements of the circuit corre- 
spond to displacement quantities. In the second method the voltages at 
nodal points correspond to displacement quantities and the currents 
correspond to force quantities. The second method will be used in all 
analogies throughout this paper since this type of analogy has been 
found to be more efficient on the Cal -Tech analog computer for beam and 
plate analogies. The types of circuits to be considered are those con- 
taining only passive elements. The elements which \d.ll be used. are 
inductors, resistors, capacitors, and transformers. (Amplifiers are not 
required in the analogies described herein. ) 

Another simple principle can now be stated for circuits with passive 
elements. In many cases it is possible to determine the correspondence 
between currents and force quantities after the correspondence between 
voltages and displacements has been assumed. Whenever the voltages at 
either end of a circuit element both represent the same type of displace- 
ment quantity, the product of current by the drop in voltage across the 
element must correspond to a true structural energy quantity (or deriva- 
tive or integral with respect to time of a structural energy quantity) . 
The structural energy quantity may be either Internal strain energy or 
external work. 

In the equations which govern the structure there are three types 
of relations Involved. Forces are related to forces from equilibrium 
conditions. Forces are related to displacements from an application of 
Hooke’s law. Displacements are related to displacements from geometric 
conditions of continuity. The equations of the electrical circuits 
involve three corresponding types of relations. Currents are related to 
currents from Klrchhoff 's nodal law which states that the outflow of 
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current at a node must equal the inflow. Currents are related to volt- 
ages from Ohm's lav for a resistor and similar laws for an inductor or 
a capacitor. Voltages are related to voltages, in the analogies herein 
described, by means of the transformer law. (The voltage relations 
given by Kirchhoff's mesh law are not employed.) The correspondence 
between these structural and" electrical relationships may be illustrated 
as given in table 1. 

A final principle may be stated in relation to the coefficients of 
the equations which govern the electrical circuit. For maximum simplicity 
in the process of solution it is necessary to design the circuit in such 
a way that the equations governing the circuit have coefficients which 
are independent of the frequency. 


EQUATIONS GOVERNING ELECTRICAL CIRCUIT ELEMENTS 


The equations governing currents and voltages in individual circuit 
elements will be stated in complete detail. These equations will be 
expressed in the form of first-order difference equations. The first 
equation to be considered will be an expression of Kirchhoff's law for 
a node. Currents at the kth node of a circuit are shovm. in figure 1(a). 
The cirrrent in the nth element of the main line is indicated as in. 

The cTurrent being forced into the line by a generator at the kth node is 
Ik. Kirchhoff's law gives 


^n “ ^-1 ” ^k ” ^ 


( 1 ) 


A differencing operator Ak will be introduced and defined as the jump 
in a function across the kth node. When this operator is applied to the 
ciirrent in the result is 


^^n = ^n - ^-1 


( 2 ) 


Substituting equation (2) into equation (l) gives 


^^n - ^k 


(3) 
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The next equation to be considered is an expression of Ohm’s law 
for a resistor as shown in figinre 1(b). This resistor, which has a 
resistance is the nth element connecting two nodal points in a 

circuit. The current in the element is related to the voltages at the 
nodes by the equation 


Rnin = % - Ek+i 


(M 


A second differencing operator ^ will be introduced and defined as 
the Jump, or increase, in a fimction across the nth element. Applying 
this operator to the voltage gives 


AnEk = Ek+l - Ek 


(5) 


Substituting equation (5) into equation (4) gives 


^nEk = -Enin 


( 6 ) 


Similar equations can be written for the inductor and condenser as shown 
in figures 1(c) and 1(d). The equation for the inductor becomes 




(7) 


The equation for the capacitor becomes 

dEv ■] 

The final element to be considered is the ideal transformer as shown 
in figure l(e). The ideal transformer is governed by two laws, one 
relating the currents in the two coils and one relating the voltage drops 
across the colls. The turns ratio r^^ of the coils, as shown in the 

figure, enters into both equations. Using a prime to indicate quantities 
in the upper coll the equations relating currents and voltages are as 
follows : 
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Ai^k* = ~ 


( 10 ) 


The currents are related by an algebraic equation -while the voltages are 
related by a difference equation, since a.n voltages are defined with 
respect to the same ground. 


BEAM WITH TORSIONAL LOADING 


A beam with torsional loading is sho-wn in figure 2(a). It -will be 
assumed that the loading consists of a finite mmber of concentrated 
torques. If the applied load is distributed, it must be replaced in 
some manner by an equivalent set of concentrated loads. An element of 
the beam at the kth station is shown in figure 2(b). The equation of 
equilibrium of this element, relating applied load and section torque 
(fig. 2(c)), is as foUo-ws: 


AkTn = -Qk 


( 11 ) 


The operator Ak indicates the junp in a function across the kth station. 

The rotation (angle of twist) at any point is related to the section 
torque by the equation 


dx GJ 


( 12 ) 


The difference between 0 at either end of the nth bay, or segment of 
the beam, is given by the area under the T/Go diagram as can be seen 
from an integration of equation (12) as follows: 



(13) 


The length of the nth bay will be indica-ted as and the average value 

of 1/GJ over the nth bay as l/GJn. Since T is constant over the 
nth bay equation ( 13 ) may be -written as follows: 
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- 0k = (5^)^n 

The 1/GJ diagram is shown, in. figure 2(d) where the average values of 
l/GJn are indicated as a stepped diagram. The rotations are shown in 
figure 2(e) as a broken-line diagram. The operator An may be intro- 
duced and defined as the increase in a function across the nth bay. 
Using this operator^ equation (l4) becomes 



The structural problem is now governed by two first-order difference 
equations given as equations (ll) and ( 15 ). It is now necessary to 
choose two electrical equations having the same. form in order to design 
the analogous circuit. Equation (ll) relates forces to forces and must 
correspond to an electrical equation relating currents to currents. The 
only electrical equation having the correct form is equation (3) which 
expresses Kirchhoff's nodal law. By comparing equation ( 3 ) with equa- 
tion ( 11 ) the correspondence between currents and forces becomes obvious. 
This correspondence is illustrated in figure 3(a). 

Equation (I 5 ) relates forces to displacements and must correspond 
to an electrical equation relating current to voltage. Equations {6), 

(7)j and (8) each have the correct form and any one of these equations 
may be used. Equation (6) involving a resistor will be chosen for the 
present case. The correspondence between electrical and structural 
quantities is shown in figure 3(b). The components of the circuit which 
are shown in figures 3(a) and 3(b) may be joined together to obtain the 
complete circuit as shown in figinre 3(c). 

In order for the analogy to be complete it is also necessary that 
the end conditions shall be correctly represented in the circuit. At 
an end -vdiere the rotation is zero the corresponding voltage is brought 
to zero by grounding the circuit* At an end where the section torque is 
zero this force condition may be taken care of directly without con- 
verting it to a condition upon the derivative of the rotation as is nec- 
essary in analytical solutions. It is necessary that the corresponding 
current shall be zero. This merely means that there must be no circuit 
elements beyond the end of the beam. Such a boundary condition is often- 
times referred to as a natural boundary condition. 

The elements of correspondence between the structure and the analogous 
circuit may be outlined as given in table 2. 
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BEAM WITH HEEDING LOADS 


The problem of torsional loading was explained in detail in order 
to provide a simple introduction to the process of designing an analogous 
circuit. The problem of bending due to vertical loads will now be 
treated more briefly. The distributed loading is assumed to have been 
replaced by a set of equivalent forces with acting at station k 

as shown in figure il-(a). Vertical equilibrium of an element at the 
kth station, as shown in figure ^(b), provides the following equation 
relating section shear to applied load: 


AkVn = -Pi 


(16) 


Equilibrium of moments acting upon the nth bay} as shown in figure ^(c), 
gives the following relation between bending moment and shear: 


^iMk = ^n^n 


(IT) 


The remaining equations may be most conveniently derived by con- 
sidering the conjugate beam loaded with the M/EI diagram. This loading 
must be assiimed to be replaced by a set of equivalent concentrations as 
shown in figure U(d). The slopes are computed as the shears in the con- 
jugate beam as shown in figure k{e). This gives the following relation 
between moment and slope: 



The quantity I/EIj^ may be computed as the average value of l/EI over 
the length which extends between midpoints of the bays which are 

adjacent to the kth station. The deflections are obtained as the bending 
moments in the conjugate beam as shown in figure 4 (f).' The relation 
between slope and deflection is given by 


- XnPii 


(19) 


The structural theory is expressed by four first-order difference 
equations given as equations (16) to (19). It is now necessary to choose 
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four corresponding circuit equations. Equation (16) relates forces to 
forces at the kth station and calls for an application of Kirchhoff's 
law at the kth node of a circuit. Corresponding currents are illustrated 
in figure 5(a) and are related hy the equation 


^k^ “ ^k 


( 20 ) 


Equation (I7) is also a relation between force quantities and calls 
for another application of Kirchhoff's law. Since equation (17) relates 
to the nth bay of the struct\n:e, the corresponding currents must be 
associated with the nth node of a circuit as shown in figure 5 (b). These 
currents will be indicated with primes and are related by the following 
equation: 


Anik' = In' 


( 21 ) 


The next structural equation to be considered is equation (I8) which 
relates slope to moment. This is a relation between a displacement quan- 
tity and a force quantity. The corresponding electrical equation may be 
obtained by considering a resistor between two nodes as shown in fig- 
ure 5 (c). Since the current in this resistor is analogous to bending 
moment the c\irrent is indicated with a prime. Ohm's law gives the fol- 
lowing equation: 


Ak^n' = -Rklk* 


( 22 ) 


One structural equation remains to be considered. Equation (19) 
relates displacements to displacements. This type of equation did not 
occur in the torsion problem. An electrical equation relating voltages 
is required. Such an equation is given by the transformer law. Comparing 
equation (lO) with equation (19) indicates that these equations are not 
of the same form. However, the electrical equation will take the same 
form as the structural equation if one coil of the transformer is grounded 
as shown in figure 5 (<i)* If the turns ratio of the transformer is rn, 
the relation between voltages is given by 


AnEk = TnEn* 


(23) 



NACA TN 2785 


11 


It is now necessary to assenible the components shown in figure 5 
into a complete circuit in which the currents and voltages will he con- 
sistent. The required arrangement is shown in figure 6. An examination 
of this figure, along with figure 5# will show that all of the voltages 
and currents are consistent. The group of four components consisting of 
two nodes, a resistor, and a transformer may be referred to as one cell 
of the circuit. Not only must it be shown that the voltages and currents 
will be consistent when the circuit elements are connected but it must 
also be shown that the current law for the transformer is not violated. 
This current law is given by equation ( 9 ). In terms of the currents 
shown in figure 6 this current law becomes 

In* = -^nin (24) 


If the corresponding structural quantities, as illustrated in figure 5; 
are substituted into equation (24), the equation becomes an Identity. 
Hence, this condition will always be satisfied. When the circuit ele- 
ments are connected the circuit appears as shown in figure 7- 

The elements of correspondence of the analogy may be siommarized as 
in table 3- 


BEANE WITH TOESIONAL AND BENDING VIBRATION 


The structural equations for a vibrating beam may be conveniently 
formulated by using D’Alembert’s principle. In the case of torsion the 
concentrated torque loafL may be expressed in terms of equivalent concen- 
trated mass moments of inertia as follows: 

The additional structural equations are the same as for the static case. 
Equation (25) expresses a relation between forces and displacements. 

When a structural equation relating forces and displacements contains 
time derivatives it is not immediately obvious as to whether the corre- 
sponding circuit element should be an Inductor, resistor, or capacitor. 
This decision can be made only after the complete circuit is designed 
and assembled. The condition of consistency of the currents, voltages, 
and their time derivatives will indicate whether or not the correct ele- 
ments have been chosen. 
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In the present case the correct element Is a condenser as shovn in 
figure 8 (a). The relation between current and voltage is given by 

Ik = -Ck (26) 

Comparing equation (26) with equation (25) it is seen that voltage is 
analogous to rotational velocity. It is possible to define an analogy 
in which voltage is analogous to rotation by merely differentiating 
equation (26) to obtain 


dt 



( 27 ) 


If equation (27) is compared with equation (25) it is seen that displace- 
ment corresponds to voltage and torque corresponds to the time derivative 
of current. Equations (26) and (27) are equally satisfactory for defining 
the analogy. Throughout this paper voltage will always be assijmed analo- 
gous to displacement and equation (27) will be used in the present case. 
The reason for this choice is that in calculating vibration inodes the 
displacements are of fundamental interest. In a transient stress prob- 
lem currents are usually made analogous to forces because the forces are 
of fundamental interest. 


The structural equation of equilibrium as given by equation (ll) 
corresponds to the following circuit equation; 



dt 


(28) 


£ 


e 


The nodal point for this equation is shown in figure 8 (b). The struc- 
tural equation relating section torque to rcjtation, as given by equa- 
tion (15)^ corresponds to the circuit element shown in figure 8 (c). In 
order that the currents and voltages in the final circuit shall be con- 
sistent, it is necessary that this element be chosen as an Inductor. The 
electrical equation is as follows: 


^n^k - “^n 


dij 

dt" 


(29) 


The complete circuit diagram is shown in figure 9 > The elements of 
correspondence of the analogy may be summarized as in table 1 . ^ 


p# 
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In the case of tending vibration the applied load is an inertia 
load given by the formula 


Pk = -“k 


dt2 


(30) 


The corresponding electrical equation is obtained from a condenser as 
shown in figure 10(a); 


dij, _ a^Ej, 

at"' 


(31) 


The remaining stnjctural equations for the vibration case axe the same 
as for the static case. 


The equation relating currents at the kth node of the shear circuit, 
as shown in figure 10(b), is as follows: 



(32) 


The equation relating ciirrents at the nth node of the moment circuit, as 
shown in figure 10(c), is as follows: 



din’ 

dt 


(33) 


The structural equation relating slope and moment corresponds to an 
electrical equation relating voltage and current in a circuit element 
as shown in figure 10(d). This element must be an Inductor in order to 
obtain consistent voltages and currents in the complete circuit. The 
electrical equation is 


AkEn’ = -Lk 


dt 




A transformer is required for the vibration case just as in the static 
case and the voltages are again related by equation (23) • 

The complete circuit is shown in figure 11. The elements of corre- 
spondence of the analogy may be summarized as in table 5« 
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BEAM IN BENDING WITH AXIAL TENSION 


A static -loading case will be considered for a beam with vertical 
loads and axial tension as shown in figure 12(a). The distributed 
loading has been replaced by a set of equivalent concentrated loads. A 
small element of the beam at the kth station is shown in figure 12(b) 
with the vertical load P^. The resultant internal forces have been 
resolved into vertical and horizontal components. The horizontal com- 
ponent is the tension Tn in the nth bay. This quantity is known or 
may be computed by direct quadrature as in the case of centrlfiigal forces. 
The vertical component Vsn. Is the statically determinant part of the 
section shear and may be obtained by direct integration of the vertical 
loads. The condition for vertical equilibrium gives 


AkVsn - -Pk 


(35) 


The total shear on a cross section consists of a statically deter- 
minate part and an indeterminate part Vj^^. The total shear may be 

expressed as the difference of the two parts 


V = 
''n 


'sn 


- V 


in 


(36) 


The total shear acts on a plane normal to the axis of the beam as shown 
in figure 12(c). Since the slope of the axis is a snlall angle it 

is seen from figure 12 that the section shear can be obtained from the 
following formula; 


Vn = Vsn - TnPn 


(37) 


Comparing equations (36) and (37) gives the formula for Vin: 


Vin = TnPn 


(38) 


The remaining structural equations are the same as for the case of 
bending without tension and are given as equations (17)^ (l8), and (19). 
It is convenient in designing the analogous circuit to substitute equa- 
tion (19) into equation (38) in order to obtain a difference equation as 
follows ; 
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= 



(39) 


The analogous circuit for the present case differs from that for 
Vxiding without tension by changes in the shear circuit. The new ele- 
ments which are required axe shown in figure 13. At the kth node the 
equation relating currents is given as before by equation (20). The 
current in, as indicated in figure 13(a), represents only the statically 
determinate part of the shear in agreement with equation ( 35 ) . Equa- 
tion (39) relates a force to a displacement and calls for a new circuit 
element in the shear circuit. A resistor must be used as shown in fig- 
ure 13(b) if a resistor is used in the moment circuit. The current 
through this element is the indeterminate part of the shear. 

Equation (36) expresses an algebraic relation between force quan- 
tities. An analogous circuit equation can be obtained by considering 
two elements in parallel as shown in figure 13 (c). One coll of the 
transformer, carrying a current which is analogous to the total shear, 
is in parallel with the resistor carrying the indeterminate part of the 
shear. The relation between currents, from Kirchhoff's law, is seen 
from the figure to be 


i ' ' ' 
■^n 


= in - in” 


(40) 


The circuit elements may be arranged as indicated in figure l4. If 
the voltages and currents which have been previously associated with each 
element are examined, they will be found to be consistent in this arrange- 
ment, The complete circuit with loading generators is shown in figure I5. 
It is interesting to note that the addition of axial tension makes only 
a minor Increase in difficulty when the solution is being obtained on an 
analog computer. However, when the solution is being obtained on a dig- 
ital computer, the addition of axial tension makes a very large increase 
in the calculation effort since the solution can no longer be obtained 
by direct quadrature. This analogy is useful in the analysis of such 
structures as helicopter blades or suspension bridges. 


EFFECT OF SHEARING STRAINS ON BENDING DEFLECTIONS 


The effect of shearing strains upon the bending deflections of beams 
may be easily included in the analogous circuit. The only structural 
equation which must be revised is that which relates slope to deflection. 
In figure l6(a) there is shown a short segment of a beam within the 
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nth. bay. A vertical line has been drawn normal to the axis of the beam 
as it would appear before loading. After loading the axis deforms into 
a curve as shown in figure l6(b). The physical points lying on the 
original normal now lie on the solid line which makes an angle with 

a perpendicular to the axis. The original normal rotates, because of 
the loading, through an angle from its original vertical position. 

Introducing 0 as the slope of the axis of the beam at any point 

gives 

0 = ^ (41) 

8x 

From figure l6(b) it is seen that 0 may be expressed as 

0 = p + 7 (42) 


It is necessary to express the shearing strain in terms of the shear by 
the formula 



(43) 


In this formula the area A is the area of a cross section which is 
resistant to vertical shear. This section property was introduced and 
defined by Kirste (reference 1) for thin-walled beams with closed cross 
sections such as airplane wings or fuselages. Substituting equation (43) 
into equation (42) gives 

9 = 3 + rr (tt) 


Introducing finite differences into equation (4l) gives 

Vn = ^^k (^5) 


Writing equation (44) for the nth bay gives 


0 


n 


= P. 


n 


GA^ 


(46) 
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Equations (4-5) and (46) are the new structural equations which must be 
represented by electrical circuits and analogous electrical equations. 

The required circuit elements to allow for the effect of shearing 
strains are shown in figure I7. The nth node in the moment circuit is 
shown in figiare 17(a). Equation (46) relates force and displacement 
quantities. The analogous electrical relations from Ohm's law for a 
resistor are illustrated in figure 17(b). The circuit equation which 
is analogous to equation (46) becomes 


En" - En' = Rn’In’ 


(i^7) 


Equation (45) relates two displacement quantities. The analogous elec- 
trical equation is obtained from the transformer law and is illustrated 
in figiire 17(c). The electrical equation becomes 


rnEn* ‘ = ^E^ 


(48) 


When the elements are assembled into a conplete circuit it is found 
that in order to obtain consistent cxurrents the current in the 

resistor shown in figure 17(b) must be analogous to the product 
rather than as would be suggested by con^aring equations (46) and 
(47). The assembled circuit is shown in figure I8. 


EFFECT OF ROTARY MERTIA UPON BENDING VIBRATIONS 


It is possible to Include the effect of rotary inertia in the cir- 
cuit for bending vibration in a rather simple manner. In order to for- 
mulate this analogy it is necessary to consider the effect of applied 
moment loads upon the behavior of a beam. In figure 19(a) there is shown 
a cantilever beam with moment loads distributed along its length. These 
moment loads are similar to the twisting moments which occur in elastic 
plates. The distributed moment load must be replaced by an equivalent 
set of concentrated couples. The nth bay of the beam is shown in fig- 
ure 19(b) with a concentrated couple T^^. Rotational equilibriiam of this 

element gives the following structural equation: 


(i^9) 


AijMk = Xn.Va + Tn 
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The remaining structural equations are not changed from those given 
previously for bending vibration. If the moment loads are due to rotary 
inertia^ D'Alembert’s principle gives the following formula for the con- 
centrated couples: 






(50) 


It is assumed that the distributed rotary inertia has been replaced by 
an equivalent set of concentrated values which are Indicated as Ip^ for 

the nth bay. 

The analogous circuit equation for equation ( 49 ) is obtained from 
Kirchhoff's nodal law. The nth node of the moment circuit is shown in 
figure 19(c). Kirchhoff's law gives 


Z^ik' = In’ - In” 


(51) 


Equation (5I) is analogous to equation ( 49 ). In order to obtain an elec- 
trical equation which is analogous to equation (50) it is necessary to 
introduce a condenser at the nth node of the moment circuit as shown in 
figure 19(d). The applied concentrated couple is analogous to the time 
derivative of the current and the appropriate electrical equation becomes 


— — — — 

dt ^ at2 - . 


(52) 


For the vibration problem it is necessary to differentiate equation ( 51 ) 
with respect to time to obtain 


^ dik’ ay ^^n" 

“ dt dt " dt 


( 53 ) 


The new elements of correspondence in the analogy become immediately 
apparent by comparing equations (52) and (53) with equations ( 49 ) and 
(50). The complete circuit is shown in figure 20 . 
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GRIDWORK OF BEAMS HAVING NO TORSIONAL STIFFNESS 


A cantilever gridwcrk; as shown in figure 21(a), will be considered. 
The structure is loaded with vertical loads normal to the plane of the 
gridwork. It will be assumed that the ribs and spars have no torsional 
stiffness. The condition of structural continuity requires that the 
ribs and spars must have the same vertical deflections at the intersec- 
tion points. There is no relation between the slopes and rotations of 
the ribs and spars. From these considerations it may be realized that 
the solution for a sweptback gridwork, as shown in figure 21(b), is 
identically the same as for the straight gridwork for a load at the same 
corresponding point. 

Since the spars and ribs act as ordinary beams in bending, they may 
be represented by analogous circuits which have been previously presented. 
One cell of the circuit for the kth rib is shown in figure 22(a) and one 
cell of the circuit for the Jth spar is shown in figure 22(b). The nota- 
tion of the spar circuit is exactly the same as was previously used for 
a beam. 

In the circuit for the kth rib the jth node in the shear circuit 
corresponds to the structural intersection point of the kth rib and 
Jth spar. Similarly, in the circuit for the Jth spar the kth node corre- 
sponds to the same structural intersection point. Consequently, from 
the structural condition of continuity the voltage Ej (corresponding to 
wj) and the voltage Ek (corresponding to wk) must be the same under 
all loading conditions. This requires that the shear circuits be con- 
nected to form a planar network as shown in figure 22(c). The volt- 
age Ejk (corresponding to wjk) is identically the same as Ej or Ek. 
Similarly, the loading current Ijk (corresponding to the concentrated 
load Pjk) is identically the same as Ij in the rib circuit or Ik in 
the spar circuit. 

The only boundary conditions idilch need to be considered in desig ning 
the circuit are the conditions that the deflection and slope of the spars 
must be zero at the root. This condition is obtained by groimding the 
shear and moment circuits at these points. The remaining boundary con- 
ditions are natiiral boundary conditions which are automatically satis- 
fied by the analogous circuit. 


ADDITIONAL ANALOGIES 


The anal ogles which have been developed were chosen as the most 
elementary types which might be used to provide an introduction to this 
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subject for the structural engineer. Many additional analogies have been 
published, although some of them do not satisfy the condition of being 
true dynamic analogies in the sense that they are valid for any time 
history of the variables. 

Analogous circuits for beams with coupled bending and torsion are 
given by McCann and MacNeal in reference 2 . This reference gives cir- 
cuits for several complete beams including bopidary conditions. Analogous 
circuits for elastic plates have been given by MacNeal (reference 3 ) • 
Applications of the beam analogies to aeroelastlc problems are described 
in reference 4 . 


SCALE FACTORS 


Once it has been determined that a mechanical system and an elec- 
trical netvork satisfy the same form of equations, the electrical net- 
work could be tonstructed by equating currents to forces, voltages to 
displacements, and electrical time to mechanical time. This procedure 
will usually lead to a group of physical electrical elements that are 
of very inconvenient magnitudes. For example, a mass of one slug is not 
at all unreasonable but the equivalent electrical element, a capacitor 
of one farad, is very unreasonable. 

To avoid such difficulties, scale factors (or modeling factors) are 
introduced in the relationships between mechanical and electrical quan- 
tities, and these factors are chosen so as to make the resulting elec- 
trical circuit as reasonable as possible as regards the availability of 
elements of the required sizes and the accuracy of measurements. 

One scale factor is introduced relating mechanical time and elec- 
trical time, and another, relating mechanical energy and electrical 
energy. A group of scale factors relating the dependent variables (dis- 
placements and voltages), in the manner of a linear transformation of 
coordinates, is also introduced. Usually the matrix of coefficients of 
this linear transformation will be a diagonal matrix so that each dis- 
placement is calculated from one voltage of the electrical circuit. In 
such cases there will be several distinct scale factors relating dis- 
placements and voltages, each of which is chosen to give the most reason- 
able values of electrical-circuit elements. In simple exanples the scale 
factor may be the same for all displacements giving a total of three 
scale factors. 

The details of this process will not be further discussed in this 
report since it relates to the design of circuits for use on analog com- 
puters and not to the form of the required circuits and the equations 
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describing them. In practice the form of the circuit Is first determined 
and values are given for all of the elements in algebraic form as if all 
the scale factors were unity. At a later stage scale factors are chosen 
and numerical values of the elements are computed. A more complete dis- 
cussion of scale factors is included in reference Jj-. 


CONCLUDING EEMARKS 


In order that electrical analogies may find their rightful amount 
of utilization in structural analysis, it is necessary that the struc- 
tural engineer be enabled to acquire a reasonable understanding of the 
process whereby an analogous circuit is designed. This will enable him 
to know where electrical analog methods are useful and what degree of 
accuracy may be expected from such solutions. An atten^jt has been made 
herein to explain in very elementary terms the procedure for designing 
circuits which are analogous to statically deteanninate beams under various 
conditions. It is hoped that structxiral engineers will be enabled to 
proceed from this paper to the reading of references which are written 
from the viewpoint of electrical engineers. En5>hasis has been placed 
upon the fact that the type of analogy described herein is a true dynamic 
analogy which permits the struct\xre to be analyzed for transient as well 
as static and steady-state conditions. 


California Institute of Technology 

Pasadena, Calif., January .9 ^ 1952 
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TABLE 1 

CORRESPONDENCE BETWEEN STRUCTURAL AND ELECTRICAL RELATIONSHIPS 


Structural 

Electrical 

Equations relating - 

Forces to forces 

Currents to currents 

Forces to displacenents 

Currents to voltages 

Displacements to displacements 

Voltages to voltages 

Equations expressing - 

Equilitrivun 

Kirchhoff's nodal law 

Hookers law 

Ohm's law 

Continuity 

Transformer law 
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TABLE 2 

ELEMENTS OF COREESPOHDENCE BETWEEN STRUCTURE AND ANALOGOUS 
CIRCUIT FOR BEAM WITH TORSIONAL LOADING 


Structural 

Electrical 

AkTn = -Qk 

Akin = Ik 


An®k = "^n.In 


Qk 

Ik 

-Tn 

in 


Rn 

QJn 

0 k 

Ek 
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TABLE k 

ELEMEMTS OF COEEESPOKDMCE BETWEEN STRUCTURE AND ANALOGOUS 
CIRCUIT FOR BEAM WITH TORSIONAL VIBRATION 


Sti^ctxiral 

Electrical 


dik ^ d^Ek 

dt ■ ^ dt2 

AfcTn = -Qk 

fdlA _ dlk 

^(dt ) dt 


V^k = ^ 

"0k 

Ck 

Ok 

dt 

0E 

Ek 

-Tn 

din 

dt 

GJn 

Ln 
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TABLE 5 -■ 

ELEMENTS OE COBBESPONDENCE BETWEEN STRUCTUBE AND ANALOGOUS 
CIRCUIT FOR BEAM WITH BENDING VIBRATION 
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(a) Nodal point. = Ifc* 


k " k-t-l 

A/WW\MAA< • 

— ; ► 

(b) Resistor. AqEjj. = -R^in. 

U U*i 

-iUinnnnnnn-4r' 

—I — » 

dij, 

(c) Inductor. = -L^ -— . 

CLu 


k 






Lh 


dEr 


(d) Condenser. An in. 


dt 


k k+! 




k-H 


(e) Transformer, i^* = -r^ij^; = — A^Ej^,. 


Figure 1.- Electrical circuit elements. 


t 
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(a) Loads. 

Qk 



(L) Beam element at kth station. 



h 


(c) Section torque. 



Figure 2 .- Beam in torsion. 
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(a) Nodal point. Ajj-Tq^ = = Ik- 



(c) Typical eectioii of circuit. 

Figure 3 .- Analogous circuit for beam in torsion. 



(a) Loads. 

r-^ 

I//,./ II II 

("b) Beam element at kth station. Ajj^V^ = 



(c) Beam element in nth hay. AnM]^ = X,. 



uxn 





(d) Loads on conjugate beam. 



J 


+ 



(f) Deflections AqWjj^ = XtiPh. 
Figure 4.- Beam in bending. 
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* A 


h-Fk 


(a) kth node in shear circuit. = -Pj^; ^k^n ” ^k* 


^kSJUk 


t 




■n ’'n-n 

(h) nth node in moment circuit. = XuV^; ^Ik' = In’ • 



/»-/ i 

a hk 


(c) kth circuit element in moment circuit. = - (W*” 

Aj?n’ = -Rk^k'' 


«• 

k - 


xJ 


'cnrrr 


+ 

+ 


k+/ 


(d) nth circuit element in sheax circ\iit. ^®k ~ ^n®n’ 

Figure 5 .- Breakdown of circuit for heam in tendi n g. 
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r T 1 




(a) Condenser. Qj, = -l0^ _j _ = 

shizl'’ 


‘’Tfr 


dt2 


(b) Nodal point. = -Q,^} ^k(^) = ^* 


dt 


-7; 


SJn 


n 


(c) Inductance 


• = W 


^n» 


Fig^^re 8 .- Analogous circuit elements for torsional vibration. 
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^k) 


k 


dt 







♦ 





Figure 9 - - Complete analogous circuit for torsional beam vibration. 
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(£/,• w^) 




(a) Condenser at kth node of shear circuit. 


Pv = — 5" 

^ ^ dt2 


d^T, dl 


.. ^ = -Cv 
" dt ^ 


d^Ej 

dt^ 



(Td) kth node in shear circuit. 


'^k^n ~ "^k» '^k 



£k 
dt ■ 



/div’\ dl ’ 

(c) nth node in moment circuit. A^I j = ^ . 


eft^M 

omnnsT'— ^ 

At 


(d) 


/Xv \ ^ik’ 

Inductance in moment circuit. Aj^P^ = ■[^^“Pk^ ^k®n' ^ ”^k ‘ 


Figure 10.- Analogous circuit elements for bending yihration. 
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(a) Loaded 'beani. 



<"/>-/) 


(b) Statically determinate shear. 



\/n-i 


(c) Total shear. 

Figure 12.- Beam in bending vith axial tension. 
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s/7 


Xl^Pu 


(a) kth node in shear circuit. = -Pj.; Aj^i^ = 

- - Vi^ 

-JL^WWWW-^^ 

fif- - 


(h) Resistance element in shear circuit. 


^ ^in^ “ ~^**^** 


— v wvvw 


ktf 


//,"=- 1/ 


//.-14 


//7 


S/7 


(c) Parallel elements in shear circuit. 


V = V 
’n ’sn 


VinJ in’” = in ' in”* 

Figure 13.- Circuit elements for heam in bending with axial tension. 
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Figure l 4 .- Exploded diagram of one cell of circuit for beam in bending 

with tension. 



o» 

TO 



Figure 15 .- Analogous circuit for beam in bending with axial tension. 
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(a) Normal to "beam aocis before distortion. 



(b) Position of normal after distortion. 
Figure l6. - Distortion of an element of a beam. 
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(a) nth node in moment circuit. = ^n^n^ ~ ^n** 

I»'-An Vn 
CEA'* e») 

(h) Resistance element required by shearing strains. 




(c) Transformer in shear circuit. = ^n®n> ^®k - ^n®n ' * * 



Figure I7. - Circuit elements to give effect of shearing strains on 

"bending deflections. 



■p- 



Figure l 8 . - Analogous circuit to give effect of ahearing straina on 

bending deflections. 
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(a) Beam, loaded with moment loads. 





(h) nth hay of beam. 


I 


In"‘-rn 


\r. 


Vl^'-P^nVn 

(c) nth node in moment circuit. + Tjj^; A^ij^* = 


I ’ - I ’ ’ 

•^n -^n • 




I 


-L. “ ^ 

dt ‘ 




(d) Condenser element for rotary inertia. 

_ T ^ c 

' 3n ^4-2 ' dt n 2 • 


Figiire 19 • - Moment loads due to rotary inertia. 
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